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Abstract 



The behavior of an ideal D-dimensional boson gas in the presence of a 
uniform gravitational field is analyzed. It is explicitly shown that, contrar- 
ily to an old standing folklore, the three-dimensional gas does not undergo 
Bose-Einstein condensation at finite temperature. On the other hand, Bose- 
Einstein condensation occurs at T ^ for D = 1,2, 3 if there is a point-like 
impurity at the bottom of the vessel containing the gas. 
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I. INTRODUCTION 



The response of quantum systems to the influence of external background fields is of 
utmost importance in a wide number of physical applications. As well, the role of disorder, 
i.e., the presence of impurities in condensed matter systems, is often crucial in the occur- 
rence of remarkable physical effects. It is the aim of the present paper to investigate the 
behavior of an ideal boson gas in the presence of a uniform (i.e., constant and homogeneous) 
gravitational field and of extremely localized (actually point-like) impurities affecting the 
quantum dynamics of the bosonic particles. 

It is well known since a long time that an ideal three-dimensional boson gas in 
free space undergoes a phase transition called Bose-Einstein condensation (BEC), in which 
a finite fraction of its constituent molecules condenses in the single-particle ground state. 
Such a condensation differs from the usual condensation of a vapor into a liquid in that 
there is no phase separation. For this reason, BEC is commonly described as a phase 
transition in momentum space — the particles condense into the |p = 0) state, which has 
a uniform spatial distribution. It is also well known that such a phase transition is no 
longer possible, for free bosons, in one and two dimensions — although in both cases it 
does occur in the presence of a point-like attractive potential H,§]]. A long standing popular 
belief [P],|j|-|T(]|| is that if the particles of a 3D ideal boson gas were placed in a (uniform) 
gravitational field, then BEC would still occur, but in the condensation region there would 
be a spatial separation of the two phases, just as in a gas-liquid condensation. 

In the present paper we study the exactly solvable quantum mechanical model of an 
ideal boson gas in D = 1, 2, 3 dimensions in the presence of a uniform gravitational field and 
of a point-like impurity formally described by a 5-function potential. In order to make the 
Hamiltonian bounded from below, so that the system may attain a state of thermodynamic 
equilibrium, we shall enclose the gas in a container with impenetrable walls. Concerning the 
mathematical description of a point-like impurity, it should be remarked that a 5-potential 
is generally ill-defined when D > 1, and some renormalization procedure is mandatory. 
Actually, the rigorous mathematical procedure to deal with point-like interactions involves 



the analysis of the self-adjoint extensions of the symmetric Hamiltonian operator [ XT | . In 
the present work, however, we prefer to follow a more informal approach [0 which is closer 
to the physical intuition, but reaches the same final result as the rigorous though more 
involved method of self-adjoint extensions |13|]. To be specific, we formally treat the contact 
interaction as a .D-dimensional (^-potential, then proceed to the renormalization procedure 
in physical terms, and finally obtain the so-called Krein's formula for the Green's function, 
from which it is possible to extract the energy spectrum of the single-particle Hamiltonian. 

In Section [TT] we prove that an ideal boson gas in the presence of a uniform gravitational 
field does not undergo BEC at finite temperature, except in the one-dimensional case. This 
implies, in particular, that in the three-dimensional case no phase separation occurs in 
the thermodynamic limit, at variance with the above quoted conventional wisdom. We also 
provide a rather general sufficient condition for the occurrence of BEC in a trapped ideal gas, 
which generalizes some results obtained by other authors |jl4|-|18|1 for power-law potentials. 
In Section |T| we show that the onset of BEC in a uniform gravitational field is made possible 
in D = 2, 3 if a point-like impurity (i.e., a 5-potential) is placed at the bottom of the vessel 
containing the gas. The reason is that the presence of the impurity entails the existence of 
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a bound state, whose energy gap with respect to the continuous spectrum is what is needed 
for the ideal gas to undergo BEC. In Section |^ we draw our conclusions, whereas some 
technical details are presented in two Appendices. 



II. D-DIMENSIONAL BOSON GAS IN A UNIFORM FIELD 

It is convenient to first analyze and discuss the impurity-free case, which turns out to 
exhibit, as we shall see below, rather surprising features. Thus, in this Section we shall 
study the quantum mechanical behavior of an ideal boson gas in the presence of a uniform 
gravitational field. The existence of a (single-particle) ground state is guaranteed by the 
presence of an impenetrable wall at the bottom of the vessel containing the gas. The single- 
particle Hamiltonian is given by 

H ( D \g) = ^ + mgx, (1) 
2m 

in which we have set 

x.= (xi,...,x D ) = (r,x), p = (p u ...,p D ) = (k,p). (2) 

The gas is supposed to be enclosed in a rectangular box of sides Li,L 2 , . . . ,L D , with its 
bottom fixed at the plane x = 0. Since we are interested in the thermodynamic limit, we 
can, without lack of generality, impose periodic boundary conditions in the x±, . . . ,xd-i 
directions and Neumann boundary condition^ at x = and x = Ld, i-e., 

if)(xi, . . . , xj + Lj, . . . , x D ) = i[)(xi, . . . , Xj, . . . , x D ), j = 1, • • • , D - 1, (3) 

d x if>(r, x = 0) = d x i/}{r, x = L D ) = 0, (4) 

and then take the limits Lj — > oo, j — 1, . . . , D. After these limits are taken, the eigenfunc- 
tions and eigenvalues of HQ D \g) read 

_ exp{(i/h)k- r} / k K\(kx + a' n ) 
^' klrj_ (27rft)( -D/2 V < Ai«) ' [b) 

k 2 

^n,k = 7^— - E g a! n , neN,ke R D -\ (6) 
2m 



where Ai(x) is the Airy function [19|, a' n are the zeros of Ai'(x), and the parameters k and 



Eg are defined as 



x The reason why we impose Neumann boundary condition, instead of the seemingly more natural 



Dirichlet one, will be explained in Section III 
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f2m 2 g\ 1/3 mg h 2 K 2 

All the zeros of Ai'(x) are negative, hence the energy levels E n ^ are positive. 

If D > 1 the spectrum is purely continuous and the corresponding improper eigenfunc- 
tions are normalized according to 

On the other hand, in the one-dimensional case the spectrum is purely discrete, the normal- 
ized eigenfunctions and eigenvalues being respectively 



rh ( x) - r^ Ai^ + Q 

E n = -E g a! n , n G N. (10) 

Let us first analyze in detail the Bose-Einstein condensation (BEC) for such a one- 
dimensional system. In the grand canonical ensemble the average number of particles N at 
temperature T and chemical potential fi reads 

oo ^ 

N = T<exp[p(E n -i,)]-V 

where, as usual, (3 = l/ksT. The criterion for the occurrence of BEC is that the average 
population of the excited states remains finite as the chemical potential approaches the 
ground state energy from below, i.e., 

lim N ex = lim — — < oo. (12) 

M T£i ^TSi ^ exp \P(E n -//)]-! 



Notice that the ground state population has been split off, that being the reason why the 
above sum begins at n = 2. The sequence of eigenvalues (0) is such that the above 
mentioned BEC criterion is satisfied. Consequently, Bose-Einstein condensation is expected 
to occur, although, in order to specify the critical temperature, it would be necessary to 
sum up the series, which, up to our knowledge, cannot be done analytically. Nonetheless, 
one can estimate the critical quantities using the asymptotic behavior of E n for large n |TP[ : 



E n = -E g a' n ~ E g [3vr(4n - 3)/8] 2/3 , n > 1. (13) 

This corresponds to a density of states of the form 

dn 1 

dE ~ 7T 



p(e)~Z~z e ; 3/2e1/2 > e » e 9 - (14) 



Since E g oc g 2 ^ 3 , as g — > the energy spectrum becomes denser and denser and the ground 
state energy approaches zero. Thus, in a weak gravitational field it is reasonable to extrap- 
olate in the continuum the density of states (|14]) down to E = 0. We can then approximate 
the series in Eq. (|T2p by an integral, and eventually obtain 
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N -~l 7 ^^-V» =4,(W ' (15) 



where At = h/ y/2Trmk bT is the thermal wavelength and g s (x) = Yl^Li n~ s x n is the Bose- 
Einstein function |T[. To obtain the critical temperature, we take the limit fi — > in Eq. 
(|l5l) and equate iV ex to the total number of particles in the gas; solving for T then yields 
the approximate critical temperature 



Tc ~ (47T) 1 / 3 — - ■ (16) 
Kb \93/2{l)J 

Below T c the fraction of particles occupying the ground state is given by 

^ =1 _^ =1 _(rf (17) 



N N \TJ 

The reasoning which led us to the conclusion that a one-dimensional ideal boson gas in 
a uniform gravitational field displays BEC can be easily generalized to higher dimensions 
and other types of potential. This is the content of the following theorem. 

Theorem 1 Suppose the single-particle energy spectrum of an ideal boson gas satisfies the 
following conditions: (i) there is a gap between the fundamental and the first excited en- 
ergy levels, i.e., Ei — Eq = A > 0; (ii) the single-particle partition function is finite, i.e., 
Z = Yl^Lodnexpi—fiEn) < oo, d n being the finite degeneracy of the n-th eigenvalue of 
the single-particle Hamiltonian. Then this gas displays Bose-Einstein condensation at finite 
temperature. 

Proof. If \x < Eq, the number of particles in the excited states is bounded from above by 

v d n exp[-/3(E n -/i)] exp(/?/i) ^ 

^ 1 - exp[-(3(E n - /i 1 - exp [-(3 (E 1 - //) ^ 



Therefore 

lim N ex < eXP( ^°] A , [Z - d eM-PE )] < oo, (19) 
m^-Eq 1 — exp(— pA) 

since, by hypothesis, Z and do are finite and A > 0. Q.E.D. 

We notice that the above statement may be generalized to some cases in which part of the 
spectrum is continuous or there are infinitely degenerate energy levels. This is done under 
the suitable introduction of the density of particles in the excited states and of the single- 
particle partition function per unit volume. Some explicit examples of this generalization 
are discussed in Ref. [|| and in Section |TTT| of the present paper. 



There are many papers that discuss the problem of Bose-Einstein condensation of an 
ideal gas confined in a power-law potential |T^-p^1 , mainly using some kind of semiclassical 
approximation. In particular, they predict that a one- dimensional gas displays BEC iff the 
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power-law potential is less confining than the parabolic one, i.e., V(x) oc x n , r\ < 2. Theorem 
1 shows that this condition is too strong: BEC occurs for any positive 77. It should be clear 
that the reason of such a discrepancy is not the semiclassical approximation per se, but the 
substitution of the discrete spectrum by a smooth density of states, which may miss some 
relevant features of the energy spectrum. 

Let us return to the problem of an ideal boson gas in a uniform gravitational field. We 
shall now consider the two- and three-dimensional cases. Due to the translation invariance 
along the transverse direction(s), the proper quantity to be discussed is the number of 
particles per unit area rS D ^ = lim^.^oo N/ L\ ■ ■ ■ L D _i. The density of particles in the excited 
states is then given by 



n 



(D) 



d D ~ l k 

D-l 



(2-nh) 



exp 



A t D E 9{D-i)/2 [exp P^Ega'j + 11 
3=1 



jj, < —E g a! x . 



(20) 



The integral in Eq. ( |2"0| ) is well defined for arbitrary D > 1 due to the condition fi < —E g a! x . 
Now, since X\m. x _,\g s {x) = 00 if s < 1, the first term of the series on the r.h.s. of Eq. ( |20|) 
diverges for D < 3 as /i — * —E g a[. Therefore, a two- or three-dimensional ideal boson gas 
in a uniform gravitational field does not display Bose-Einstein condensation at T 7^ 0. 
Some remarks are in order here: 

(a) At first sight, Eq. (p0|) seems to imply absence of BEC in D = 1 too. It should 
be noted, however, that in one dimension there is no integration over transverse momenta. 
Hence, in order to remove the contribution of the ground state from the sum over states in 
Eq. (pUD, one has to begin it at j = 2. Then nix (= iV ex ) has a finite limit as jj, —y —E g a! x . 

(b) It is easy to see that the absence of BEC in a two- or three-dimensional ideal boson 
gas in a uniform gravitational field in the x-direction is due to the quantization of the 
motion in that direction. Thus, any potential V that depends only on x, and such that the 
one-dimensional Hamiltonian 



2m 



+ V{x) 



(21) 



has a discrete spectrum, will do the job of hindering BEC in D = 2, 3. 

(c) There are claims in the Literature [1|,^-TD] that a three-dimensional ideal boson gas 
in a uniform field may undergo BEC at T 7^ 0. This is an artifact of approximating the sum 
in Eq. (POD by an integral (remember that Eq. ( |2"0"| ) holds true for D > 1). Indeed, using the 
density of states given by Eq. (|14|) we obtain 



E 9(D-m [exp P{E g a' 3 +/*)]«- E g ^ / dE E^ £ 



i=i 



n=l 



n 



(D-l)/2 



-(/?£,r 3 / 2 r(3/2)E 

7T z — ' 



e n(3ii 



(D+2)/2 



n=l 



n 



(22) 
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Inserting this result into Eq. (p0[), one would be led to the incorrect conclusion that BEC 
occurs at finite temperature in D = 2 and D = 3 in the presence of a uniform field, because 
lim^o 9(D+2)/2(e (Sfl ) < oo if D > 0. 

(d) It should be clear by now that none of our conclusions so far depends crucially on 
the use of Neumann boundary condition. They would remain correct, at least qualitatively, 
had we used Dirichlet or Robin boundary condition instead. 



III. D-DIMENSIONAL BOSON GAS INTERACTING WITH A POINT-LIKE 
IMPURITY AT THE BOTTOM OF THE CONTAINER 



In this Section we finally come to the most interesting physical case in which, in addition 
to the gravitational field, there is a point-like impurity at the bottom of the vessel containing 
the gas. As we shall show here, such an impurity is enough to restore BEC in the three- 
dimensional case — and to allow its existence in the two-dimensional case, in which it is 
absent with or without the gravitational field. The single-particle Hamiltonian takes now 
the form 

H {D \g, \ D ) = ^- + mgx + \ D <^(x) = H { D \g) + A D 5^(x). (23) 

Our main task will be to show that the (^-potential creates a bound state in the two- and 
three-dimensional cases, thus paving the way for the occurrence of Bose-Einstein condensa- 
tion, at variance with the impurity-free situation discussed in the previous Section. 
Our basic tool to tackle this problem is the Green's function 

(z; x, x') = (x | [H^ (g, X D ) - z] | x'} , z G C, (24) 

from which it is possible to extract the energy spectrum. A formal expression for 
G^ D \z; x, x') can be obtained by solving the Lippmann-Schwinger integral equation, 

(z; x, x') = G { D) (z; x, x') - J ' d D y g[ d) (z; x, y ) V(y) (z; y, x') , (25) 

where G { D) and G^ are the Green's functions associated to H ( Q D) and = H ( D) +V(x), 
respectively. For l^(x) = Ad ^^(x) the integral in Eq. (^) can be done trivially, resulting 
in 

G {D) {z- x, x') = G[ D) (z; x, x') - X D G { D) (z; x, 0) G {D) (z; 0, x') . (26) 

If we now set x = 0, we obtain an algebraic equation for G^ (z; 0, x'). Solving that equation 
and inserting the result into Eq. (PB|), we end up with 



G^ ( z; x , x') = Gf > (z; x, x') - G o D \^^°)G { D \z;0,x') (2?) 

As we shall see below, Gq D \z; 0, 0) is formally divergent for D > 2, but one can still give a 
well defined meaning to Eq. (|2~7|) by renormalizing the coupling parameter Ad- The resulting 
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expression, which then makes sense also for D = 2, 3, is known as the Krein's formula |ll]] an d 
encodes the one-parameter family of self-adjoint extensions of the symmetric Hamiltonian 
operator HQ D \g). This precisely corresponds to the mathematically rigorous description of 
the 5-potential. 

To complete the construction of we still have to obtain the Green's function in the 
absence of the impurity. This is done in Appendix [A], with the result 



7TK 



d 



D-l, 



(2nh) 



D 



- exp 



h 



M [£Q£<)M£(^>)] 

Ai'[e(o)] 



where the functions u(£) and v(£) are defined in Eq. ( |A7|) , £(x) is defined in Eq. 
x <{ x >) — min(max){x, x'}. Setting x = x' = in Eq. (|28| ) we formally obtain 



Gf>(*;0,0) = - 



d^k Ai [(k 2 /2mE g ) - (z/E g 



E g J (2irH)°-i Ai' [(kV2m^) - (z/E g )] 



-C 



D 



dyy 



(d-3)/ 2 A Kv - 

Ax'(y - C) 



where 



C 



D 



k u (4tt 



,(l-0)/2 



E g Y[{D -l)/2\ 



E n 



(28) 
41), and 



(29) 



(30) 



It follows from the asymptotic behavior of the Airy function Ai(a;) for large x [19|, 



M(x) 



x—>oc 



2 v ^Fx 1 /4 



exp 



x 3 / 2 [l + 0(a;- 3 / 2 )] 



(31) 



that the integral in Eq. (|29"|) diverges in the UV region for D > 2, as anticipated. (The 
integral is finite in the IR for D > 1.) 

Before we show how to make sense of Eq. fl27D for D = 2, 3, let us discuss the one- 
dimensional case, which does not need renormalization. In this case, the energy spectrum 
can be obtained by solvingF] 



- + G«(z;0,0) 



0. 



(32) 



or, more explicitly (see Appendix |A|), 



1 k Ai {-z/E g 



A X Eg Ai' (-z/Eg) 



(33) 



2 One can easily check that the residue of Gq (z; x, x') at z = —E g a' n cancels against the residue of 
the second term on the r.h.s. of Eq. at the same pole. Therefore, all the poles of G^\z; x, x') 
are given by the solutions to Eq. ( p^ ) . 
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This equation is equivalent to the imposition of Robin boundary condition at the origin, i.e., 
^'(0) + c0(O) = 0. It interpolates between the Neumann boundary condition, for Ai — > 0, 
and the Dirichlet one, for Ai — > oo. Any of these boundary conditions prevents the flow 
of particles across the origin, so any of them can be used to represent an impenetrable 
wall at the bottom of the container. Nevertheless, it is more convenient to impose Neumann 
boundary condition in the impurity-free case, because it is then possible to model an impurity 
at the bottom of the container by a 5-potential. This would not be possible had we imposed 
Dirichlet boundary condition instead. In any case, the energy spectrum obtained by solving 
Eq. fl3"3"| ) will be purely discrete and bounded from below. As a consequence, we can say that 
in the one-dimensional case the Bose-Einstein condensation actually occurs at the lowest 
discrete energy level, although the ground state energy itself as well as the critical quantities 
are shifted with respect to the previously discussed impurity- free case. 

Let us now discuss the two- and three-dimensional cases. In order to make sense of 
the denominator in Eq. fl27|), we first have to regularize G D \z; 0, 0). We shall do this by 
introducing a UV cutoff in Eq. fl2"9|), namely, 



G { D \z;0,0)^G { Q D \A,z;0,0) = 
We now add to Gg (A, z; 0, 0) the integral 



-C D I dyy 



(D-3)/2 



Ai(y - C) 
Ai'Ciz-C) 



(34) 



[ D (\..z.n) = -C v l" ' dyy {D -^ 2 (y + a )- 1/2 
'o 



a > 0. 



(35) 



It follows from Eq. ([H]) that 



Ai(y-C) 
M'(y - C) 



y^oo 



-(y-Cr 1/2 + 0[(y-C)- 2 ] 
-y 1 2 + O (Cy- 3/2 ) ; 



(36) 



hence, the integrand of Gq(A, z; 0, 0) + Jd(A, z, a) behaves like y^ D 6 ^ 2 for large y. This 
allows us to remove the UV regulator (i.e., to take the limit A — > oo) for D < 4. At the same 
time, since we have added Id to Gq D \ we must subtract it from A^ 1 in order to keep the 
combination A^ 1 + G Q (z\ 0, 0) unaltered. We may then define the renormalized coupling 
parameter A^ as 



lim 

A— >oo 



A 



- I D (A,z,a) 



D 



(37) 



where it is understood that A^ depends on A in such a way that the limit exists. We then 
finally arrive at a meaningful expression for the Green's function G^ D \z; x, x') for D = 2, 3, 
in which the denominator of Eq. ( P7| ) is replaced by the finite expression 



1 



A D 



dyy 



(D-3)/2 



A% - C) 
Ai'(y - C) 



a 



-1/2 



(38) 
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It is possible to show (see Appendix ||) that, for any finite value of A^, g D ((, a, A^) 
has a single zero Co i n the interval — oo < Co < —a[. In physical terms, this means the 
existence of a bound state with energy E = E g ( . The rest of the energy spectrum forms 
a continuum starting at E = —E g a! x . The presence of this gap in the energy spectrum is 
enough to guarantee the occurrence of BEC. The proof of this fact is similar to that of 
Theorem 1, the only difference being that what saturates in the limit [i — > E is not iV ex , 
but ni^ . Some examples of this phenomenon are discussed in detail in Ref. ||, where it is 
also shown how to obtain the critical quantities. Working in close analogy, one can obtain 
an estimate of the critical quantities in the present situation, taking Eq. (P2|) suitably into 
account. If the energy gap created by the impurity is much greater than the energy splitting 
due to the gravitational field, i.e., A = —E g a[ — E ^> —E g a' 2 + E g a[, one can obtain a good 
approximation to the critical temperature T c by solving the equation 



a£-V d > 



4tt (k\ Tc ) 3 9{d+2)/2 [exp(-A/fc B T c )] . 



(39) 



It is worthwhile to stress that now, because the bound state energy Eq is strictly below the 
continuum threshold (—Ega^), we can safely use Eq. fl22|) to estimate the critical quantities 
in£> = 2,3. 

We close this section with a somewhat technical remark. Aside from being positive, 
the parameter a in Eq. (|38D is arbitrary, and has to be fixed by some renormalization 
prescription. One possibility is the so called Bergmann-Manuel-Tarrach p0[ renormalization 



prescription, in which the bound state energy Eq labels the one-parameter family of self- 
adjoint extensions of the symmetric Hamiltonian HQ D \g). Then Eq. (EM) becomes equivalent 
to the pair of equations 



Sd(C,Co)|bmt = C d dyy 



,(D-3)/2 



Ai(y-Co) Ai(y-C) 
Ai'(y-Co) Ai'(y-C) 



(40) 



a 



C 



D 



dyy 



(D-3)/2 



Ai(y - Co 
Ai'(y - Co 



+ (y + a) 



-1/2 



(41) 



where Co = E /E g < —a[. The parameter a > is thus the subtraction point at which the 
"running" coupling parameter A^ is defined. 



IV. CONCLUSIONS 

In this paper we have explicitly solved the quantum dynamics and studied the ther- 
modynamic equilibrium of an ideal D-dimensional boson gas in the presence of a uniform 
gravitational field and a point-like impurity at the bottom of the vessel containing the gas. 
For convenience, in the present analysis we have imposed Neumann boundary condition at 
the bottom of the container, but our results might be generalized to Dirichlet or Robin 
boundary conditions without any substantial modification in the physical behavior. In the 
impurity-free case it has been shown that Bose-Einstein condensation at finite temperature 
is possible only in the one- dimensional case and an estimate of the critical temperature in 
this case has been obtained. It has also been elucidated why the conventional wisdom that 
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BEC (with a phase separation) might occur in the three-dimensional case does actually fail: 
the reason eventually lies in the illegitimate use of a continuous approximation to the density 
of states in the computation of the average number of particles in the excited states. 

On the other hand, it has been proved that the presence of a point-like impurity is enough 
to allow BEC at T ^ also in two and three dimensions. The reason is that the impurity 
creates a bound state in the single-particle spectrum, where particles can now accumulate. 
It should also be emphasized that a ^-potential in the presence of a uniform field is always 
attractive in two and three dimensions, irrespective of the sign of the renormalized coupling 
parameter. 

The main interest in the study of the present model is in its exact solvability. Nonetheless, 
it is evident that the key physical features here exhibited will persist even if more realistic 
impurity potentials are used. The situation is less clear if one considers an interacting boson 
gas (for the general definition of BEC, applicable to this case, see Ref. [f2lH ). It is reasonable 
to assume that our results still hold if the mean field interaction between the particles in 
the gas is smaller than (i) the energy splitting due to the gravitational field, and (ii) the 
energy gap created by the impurity (if the latter is present). This condition, however, is 
likely to be violated as more and more particles accumulate in the lowest energy level, until 
the interaction between the particles cannot be neglected anymore. What happens then 
awaits further investigation. 
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APPENDIX A: 

The Green's function Gq D ^ (z; x, x') satisfies the partial differential equation 

Gi D) (^x,x') = 5 (D) (x-x). (Al) 

We can reduce Eq. ( |A1| ) to an ordinary differential equation by Fourier transforming in the 
transverse coordinates: 

+ ^ + mgx ~ z ) G{ - Z) k; X) ^ = 6<yX ~ (A2) 

the Green's function will then be given byf] 

Gf\z\x.,y{) = j jT^j^ZT expj^k- (r - r')\ g(z,k;x,x'). (A3) 



3 In the one-dimensional case we have instead G (z; x, x') = Q(z, k = 0; x, x'). 
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Upon the change of variable 

P = k.t. + Fr l 

\2m 



^ = KX + E- l (^-z) ) (A4) 



Eq. ( |A2j ) becomes 

d 2 \ . , /v 



(H ,-^G^O = -^-0- (A5) 

When £ ^ Eq. ( |A5| ) reduces to the Airy differential equation. Its solution must satisfy 
Neumann boundary condition at x = 0, i.e., d^Q{^, £')U=o = 0, and it must vanish at infinity 
lim^oo = 0. Thus, 

o = Cx «(0 0(C - o + c 2 «(0 0(£ - 0. (A6) 

where 6*(x) is the Heaviside step function and 

u{i) = Bi'(£ ) Ai(0 - Ai'(&) Bi(0, v(0 = Ai(0, (A7) 
with £ = £(x = 0). To fix the constants C\ and C 2 , one imposes continuity of (?(£,£') at 

e = e, 

^' + 0,0 = ^-0,0: (A8) 
and a jump in c^£? (£,£') at the same point, 

<W + 0, - <W - 0, = — , (A9) 

^9 



obtained by integrating Eq. ( |A5| ) from £' — e to £' + e and letting e J, 0. Applying conditions 
( |A8| ) and (|A9|) to the solution (|A6|), and using the fact that the Wronskian of and 
is given by 



W{u(0, = -Ai'(eo) W{Bi(0, Ai(0} = - Ai'fo), (A10) 

7T 



we finally obtain 



— E ff Ai'(6) ' (AH) 



where £<(£>) = min (max){£, £'}. Substituting ( [All]) into Eq. ( |A3| ) gives us the desired 
integral representation of Eq. (|28|) for G^\z] x, x'). 
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APPENDIX B: 



Here we show that g D {C, A^) nas one ( an d only one) zero in the interval — oo < C < 
—a[. Indeed, for ( large and negative we may use the first line of Eq. fl36| ) to evaluate the 
integral in Eq. (|38D , obtaining 

g 2 (C, a, A*) p - C 2 In f-i") , (Bl) 



2 



a 



A? 



-2C, 



C 



fB2) 



In both cases, lim^_oo g D (^, a, A^) = — oo. On the other hand, the integral in Eq. 
becomes divergent at the origin for D < 3 if £ | ~~ a 'ii as 



Ai'(y + a[ 



Ai(ai) 
Ai"(ai) y 



<2/ 



(B3) 



(The last equality is a consequence of Airy differential equation.) Since a[ < 0, it follows that 
lim^|_ a ' g D ((,ct, A^) = +oo (D = 2,3). By continuity, we may conclude that g D ((,a, \p) 
vanishes at least once in the interval — oo < C < To show that it vanishes only once, it 
suffices to prove that g£)(C a, A§) is a monotonically increasing function of ( in that interval. 
This follows from the identity 



-^g D ((,a,\%) = E g — 
= E n (0 



1 



— + G { n D \z;0,0 

„ (g)-z\ 2 o 



H, 



(B4) 



It shows that d^g D ((, a, A^) > if z is real and does not belong to the spectrum of HQ D \g). 
This occurs, as we have seen in Section |I[ for z < —Ega[, or ( < —a[. 
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